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Abstract — We present the results of multi-frequency polarization VLA observations of radio sources from the 
complete sample of northern, radio-bright BL Lac objects compiled by H. Kiihr and G. Schmidt. These were 
used to determine the integrated rotation measures of 18 sources, 15 of which had never been measured previ- 
ously, which hindered analysis of the intrinsic polarization properties of objects in the complete sample. These 
measurements make it possible to correct the observed orientations of the linear polarizations of these sources 
for the effect of Faraday rotation. The most probable origin for Faraday rotation in these objects is the Galactic 
interstellar medium. The results presented complete measurements of the integrated rotation measures for all 
34 sources in the complete sample of BL Lac objects. 
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1. INTRODUCTION 



>q BL Lac objects are an extreme class of active galac- 
^ tic nuclei characterized by a number of peculiar proper- 
ties: weak or absent optical line emission, rapid and 
[- — strong variability in all parts of the electromagnetic 
Q spectrum, a high degree of linear polarization, and 
CO compact radio and pointlike optical structure. 

In 1990, Kiihr and Schmidt [1] compiled a complete 
^ sample of radio-bright BL Lac objects in the Northern sky, 
Q_i consisting of 34 sources satisfying the following selection 
■ criteria: (1) total flux at 5 GHz / > 1 Jy at the epoch at 
£~ which the sample was compiled, (2) two-point spectral 
^ index between 1 1 and 6 cm a! l-6 > -0.5 (5 ~ V + "), and 
C3 (3) equivalent width of emission lines W < 5 A. The 
sample included sources brighter than 20 m in the optical 
. — , with declinations 8 > 20° located outside the Galactic 
^ plane (\b \ > 10°). We are currently analyzing multi-fre- 
H quency, polarization-sensitive Very Long Baseline Inter- 
ns ferometry (VLBI) data for all sources in the sample. 

VLBI polarization observations can yield radio 
images of compact sources such as BL Lac objects with 
angular resolution of the order of a millisecond of arc. 
Such observations contain information about all four 
Stokes parameters, i.e., both the total intensity and the 
linear and circular polarization of the radiation. The 
emission of BL Lac objects is characteristically linearly 
polarized, since it is synchrotron radiation. Thus, polar- 
ization-sensitive observations carry information not 
only about the energetics and kinematics of the sources 
observed, but also about the degree of order and orien- 
tation of the magnetic fields in their jets, which are very 
important when analyzing the physical conditions in 
these objects. 

There are a number of factors that can complicate 
deduction of the direction of the magnetic field in a par- 
ticular region in a source. Three properties of the region 



must be known: first, the orientation of the electric vec- 
tor of its linear polarization, second its optical depth, 
and finally, its rotation measure. 

The first two properties flow directly from any 
multi-frequency polarization observations, and yield 
the orientation of the magnetic field in the radiating 
region relative to the observed polarization electric vec- 
tor. When a sufficiently high degree of polarization is 
detected at even one frequency, this indicates that the 
radiating plasma must be optically thin, since optically 
thick radiation cannot be more than =10-15% polar- 
ized [2]. If the plasma is optically thin, the magnetic 
and electric fields are mutually perpendicular, whereas 
the two fields are parallel if the plasma is optically 
thick. The core of the source is usually considered to be 
optically thick, while the jet becomes increasingly 
transparent — optically thin — with distance from the 
core. This is demonstrated by the presence of highly 
polarized jet components [3] (isolated regions of 
brightening of the total intensity), whose degree of 
polarization can reach 60-70% — the limit for optically 
thin synchrotron radiation [2]. The high degrees of 
polarization that are sometimes observed in the VLBI 
jets of BL Lac objects testify that the magnetic fields 
are well ordered in these regions. As a rule, the jet emis- 
sion has a steep spectrum, confirming that we are deal- 
ing here with an optically thin plasma. 

Faraday rotation of the plane of linear polarization 
occurs when polarized radiation passes through a 
region containing thermal plasma and a magnetic field 
with a non-zero component along the line of sight. The 
rotation measure can only be derived from multi-fre- 
quency polarization observations. However, wide-band 
polarimetry (1.4, 5, 8, 15, 22 GHz, for example) is not 
suited to this purpose. As is discussed in [5], the best suited 
observations are those made with the NRAO VLA at 
1.3-1.7 GHz. When planning observations designed to 
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Fig. 1. Total-intensity maps of 01 19 + 1 15 at 1.41 GHz (left) and 1.63 GHz (right) with linear polarization vectors superposed. The 
peak flux densities are 892 and 891 mjy, respectively, and the lowest contours are 0.35% left and 0.25% right of the peaks, respec- 
tively. The horizontal and vertical scales are in arcseconds. 
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Fig. 2. Total-intensity maps of 0235 + 164 at 1.41 GHz (left) and 1.66 GHz (right) with linear polarization vectors superposed. The 
peak flux densities are 833 and 1032 mjy, respectively, and the lowest contours are 0. 125% of the peaks. The horizontal and vertical 
scales are in arcseconds. 



derive the rotation measures of compact extragalactic 
sources, several requirements on the frequency interval 
to be used should be taken into consideration. First, 
there should be a large enough spacing between the cor- 
responding wavelengths (more precisely, the squares of 
the wavelengths) to detect the rotation of the plane of 
linear polarization. Second, these wavelengths should 
be long enough not to give rise to large ambiguities 
when calculating the rotation measures. Finally, the 
wavelengths should not be very different from each 
other, since otherwise, there is the possibility of detect- 



ing radiation from electrons with qualitatively different 
energies — i.e., from physical different regions in the 
source — at the different wavelengths. 

2. OBSERVATIONS 

We obtained VLA observations (C or D configura- 
tion) during four sessions: January 19, 1996 (Session 1), 
January 29, 1996 (Session 2), February 7, 1996 (Ses- 
sion 3), and June 26, 1998 (Session 4). The 1998 obser- 
vations were conducted in order to refine our initial 
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Fig. 3. Coverage of the u-v plane for the VLA observations of 0119 + 115 (left, three scans) and 0235 + 164 (right, six scans). The 
horizontal and vertical scales are thousands of wavelengths. 
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Fig. 4. Dependence of the polarization position angle % (in degrees) on the square of the wavelength (in cm ) for the sample sources. 
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determination of the rotation measure for 0119+ 115, 
since the 1996 data for this source were of compara- 
tively low quality, leading to a large uncertainty in its 
rotation measure. All the observations were carried out 
at six frequencies (1664.9, 1635.1, 1515.9, 1467.7, 
1417.5, and 1365.1 MHz), with bandwidths of 50 MHz 
for the observations at 1664.9 and 1635.1 MHz and 
25 MHz for the remaining frequencies. The data were 
calibrated in the AIPS (Astronomical Image Processing 
System, NRAO) package. The primary flux calibrators 
were 3C 286 (Sessions 2 and 3) and 3C 48 (Sessions 1 
and 4). The polarization position angles were calibrated 
using 3C 286 (Sessions 2 and 3) and 3C 138 (Sessions 1 
and 4), for which we assumed the polarization position 
angles % = + 33° and % = -9°, respectively, at all fre- 
quencies, since these sources do not have significant 
rotation measures [8]. The instrumental polarizations 



were determined using sources with observations span- 
ning a large interval (>90°) in parallactic angle — 
2150 + 173 (Session 1), 0917 + 624 (Session 2), 
1147 + 245 (Session 3), and 0235 + 495 (Session 4). 
The unpolarized sources 2352 + 495 (Session 1, m L < 1 %), 
1323 + 321 (Sessions 2 and 3, m L < 0.6%), and 0134 + 329 
(Session 4, m L < 0.8%) were used to check the results 
of the solutions for the instrumental polarizations. 

3. RESULTS 

Structure on VLA (kiloparsec) scales was detected 
in the two sources 0119 + 115 and 0235 + 164. Figures 1 
and 2 show images of the distributions of total intensity 
for these sources, with linear-polarization vectors 
superposed. The remaining sources were unresolved at 
the sensitivity level reached in our observations, and so 



70 c 



60 



50 



0745 + 245 J - 

L-r r "' 



X 

0° 



-50 



-100 k 



0820 + 225 



A 



,1 




1147-245 J. 

A'' 

-'f 



.1- 



0° 



30 



20 



1334- 127 



o° 

-20 
-40 
-60 



1732 + 389 



300 



350 



400 



450 



500 300 

A 2 , cm 2 



350 



400 



450 



500 



Fig. 4. (Contd.) 
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their images are not shown. The images in Figs. 1 and 
2 were made at different frequencies, and variations in 
the polarization angles are clearly visible. The orienta- 
tion of these vectors corresponds to that of the polariza- 
tion electric vectors, and their length is proportional to 
the degree of polarization. The sizes and shapes of the 
antenna beams are indicated in the lower left corner of 
each image. Figure 3 shows the u-v coverages for 
0119+1 15 and 0235 + 164 for our observations, which 
are typical for VLA "snap-shot" observations of 
sources at these declinations. 

Figure 4 presents the dependence of the polarization 
position angle % on the square of the wavelength for the 
1 8 sources. In all cases, these plots show a linear depen- 



dence on X 2 . The desired rotation measure is the slope 
of the line 

X (l 2 ) = Xo + RMX 2 , RM = 0M^N e BcosQdl, (1) 

where N e is the density of thermal electrons (cm 4 ), 
B the magnetic field (iiG), the angle between the direc- 
tion of propagation of the radiation and the magnetic-field 
direction, and dl an element of path length from the source 
to the observer (pc). This expression for the rotation mea- 
sure is valid for quasi-longitudinal propagation (so that the 
angle 9 is small) and a weak magnetic field [4]. 

As a rule, errors in % are several degrees, and are 
presented at the 2a level; we had several observations 
for each source at each frequency, enabling us to estimate 
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Integrated rotation measures for all sources in the complete sample of BL Lac objects 
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For a uniformly expanding Universe (q Q = 0.5) with H Q = 100 km s Mpc 
Author does not present errors for this measurement. 
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Fig. 5. Dependence of the absolute value of the rotation measure on the absolute value of the Galactic latitude. For sources with 
several RM measurements, we chose the RM value obtained at the latest epoch. The new value for the Galactic contribution to the 
rotation measure of BL Lac found by Reynolds et al. [9] is shown by the hollow circle. 
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Fig. 6. Dependence of the degree of polarization on wavelength using six sources as examples. There is no evidence for depolarization. 
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Fig. 7. Dependence of the relative degree of polarization 
m/nii at various wavelengths on the rotation measure, where 
m ; is the intrinsic degree of polarization in the absence of 
internal Faraday rotation and m is the degree of polarization 
expected for a plane-layer model [9] in which a substantial 
fraction of the rotation measure is internal to the source. The 
dashed curve shows the theoretical dependence of the ratio 
of the degrees of polarization /m^ for wavelengths Xk 
18 and 22 cm predicted by this model. 



the errors in % by comparing these independent mea- 
surements. We determined linear least-squares fits to 
the observed % dependences. Note that the polarization 
angle is determined with accuracy to within ±n% (n = 
0, 1, 2, ...), which can lead to ambiguity in the rotation 
measure RM, since we are dealing with observations at a 
discrete set of frequencies. The square of the minimum 
difference between neighboring wavelengths will deter- 
mine this ambiguity, which in our case is ±1190 rad/m 2 . 
Nevertheless, the presence of six observing frequencies 
and their proximity to one another enables us to determine 
the rotation measures with certainty. 

The integrated rotation measures for the 18 sources 
from the complete sample of Kuhr and Schmidt [1] are 
presented in the table, together with results for the 
remaining 16 sample sources obtained in earlier obser- 
vations [5-7]. These values are primarily in the interval 
|RM| = 0-100 rad/m 2 , which coincide with the Galactic 
rotation measures derived from observations of 555 radio 
sources [8]. Only two objects have |RM| > 200 rad/m 2 . It 
is probable that the Faraday rotation we have detected 
occurs primarily when the radiation passes through the 
interstellar medium of the Galaxy. Further evidence for 
this is provided by the distribution of the |RM| values for 
all the sample sources relative to the absolute value of 
the Galactic latitude, where we can see a tendency for 
|RM| to decrease with increase \b \ (Fig. 5). 

Figure 5 does not show the two highest values of |RM|, 
for the sources 0119 + 115 and 2200 + 420 (BL Lac). 
Results obtained by Reynolds et al. [9] (personal com- 
munication) based on multi-frequency VLBI observations 
show that the Galactic contribution to the observed inte- 



grated rotation measure of BL Lac is -1 10 ± 19 rad/m 2 , 
indicated in Fig. 5 by the hollow circle. We can see that 
this new value agrees well with the general tendency 
shown by the remaining rotation measures of the sample 
sources. It is not ruled out that the rotation measure for 
0119+115 also includes an appreciable contribution from 
regions near the source. 

In addition, we do not detect evidence for depolar- 
ization (Fig. 6); the degree of polarization remains essen- 
tially constant with frequency. Depolarization arises when 
the thermal plasma associated with the RM is mixed in 
with the source [10]. In this case, the polarizations from 
the far and near regions of the source add incoherently, 
and the resulting observed polarized decreases as the 
Faraday depth increases. This is shown for a plane- 
layer model [10] in Fig. 7, which presents theoretical 
depolarization curves for various wavelengths when an 
appreciable fraction of the rotation measure is internal 
to the source. 

Moreover, another piece of evidence supporting a 
Galactic origin for most of the observed rotation mea- 
sures is that the RM values apparently do not vary with 
time, as indicated by the values in the table for sources 
for which there are observations at more than one 
epoch; note that all the results in the table were 
obtained using the same observing and reduction tech- 
niques. In addition, as was shown in [11], the rotation 
measures remain constant over periods of decades, in 
spite of the presence of appreciable variations of the 
source polarizations. 

4. CONCLUSIONS 

We have derived integrated rotation measures for 
18 sources, completing such RM measurements for all 
34 BL Lac objects in the complete sample, using simul- 
taneous VLA polarization observations at six frequen- 
cies below 2 GHz. The mean value of |RM| for the 
entire sample is -40 rad/m 2 . The observed RM values 
are associated with thermal plasma external to the 
source, and probably arise when the radiation passes 
through the interstellar medium of the Galaxy. 

However, such overall modest values for the inte- 
grated rotation measures do not rule out the presence of 
substantially higher RMs and sharp RM gradients in the 
inner regions of the compact VLBI jets of a few of these 
sources. Our multi-frequency VLBI polarization observa- 
tions, currently being analyzed, will enable us to iden- 
tify such cases of non-uniform Faraday rotation in the 
immediate vicinities of the compact radio sources. 
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